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ABSTRACT 

We report the discovery of a third planetary mass companion to the GO 
star HD 74156. High precision radial velocity measurements made with the 
^ ■ Hobby-Eberly Telescope aided the detection of this object. The best fit triple 

Keplerian model to all the available velocity data yields an orbital period of 347 



O 

o 

(N 



days and minimum mass of 0.4 Mjup for the new planet. We determine revised 
orbital periods of 51.7 and 2477 days, and minimum masses of 1.9 and 8.0 Mj^p 
• respectively for the previously known planets. Preliminary calculations indicate 

^ ■ that the derived orbits are stable, although all three planets have significant 

^ ■ orbital eccentricities (e = 0.64, 0.43, and 0.25). With our detection, HD 74156 

^ ' becomes the eighth normal star known to host three or more planets. Further 

study of this system's dynamical characteristics will likely give important insight 
■ to planet formation and evolutionary processes. 

i/-^ I Subject headings: planetary systems - stars: individual (HD 74156) - techniques: 

^ ! radial velocities 

o 

O : 1- INTRODUCTION 



X 



More than 200 planets have been discovered around nearby stars other than the Sun us- 
g the radial velocity methocj^. Despite the high number of known exoplanets, the hunt for 
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Encyclopaedia website http : //exoplanet . eu 
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additional planets around nearby stars with the radial velocity method remains a fundamen- 
tal research area. Planet candidates around nearby stars are the most s uitable for foUowup 



study like photometric monitoring for transits fe.g. iHenry et al. 



200' 



Gillon et a 



McArthur et al. 



200 



Deming et a. 



ik 



200G 



Charbonneau et al. 



2000l : ICastellanol[200ol : |Bouchv et allboosl : ISato et al.ll2005l : lshankland et al.ll20od:lL6pez-Morales et al. 



2007), astrometric pertur b ation measuremen ts (e.g. [Benedict et al.ll2002l : 



Benedict et al. 



iooj [Harrington et al.l 



2006 



2006: Bean et al. 2007). thermal emission searches 



Richardson et al.l 



CowaneLalJl2Q07|), atter npted direct imaging detection (e.g. 



ical characterization (e . g.lLaughlin &: Adam s 1999: IChiang &: Murray 



2007 



Janson et al. 



Knutson et al.ll2007 
2007h. and dynam- 



20021 : lAdams fc Laughhn 



20061 : iFord et al.ll2005l : iBarnes fc GreenbergI 120061 ). The results of these types of investiga- 
tions improve our understanding of planet formation and evolution. 

Radial velocity planet searches are also a crucial component of the quest to determine 
the uniqueness of our own solar system. It should be noted that only recently have some 
surveys achieved the necessary precision and timespan to detect true Jupiter analogs (i.e. 
planets with roughly the same mass and orbital semimajor axis), and it is not yet possible to 
detect analogs of any of the other planets in our solar system. Therefore, the continuation 
and improvement of radial velocity planet searches should be a top priority. 



After our successful detection of an additional companion to Cnc (= 55 Cnc. lMcArthur et al. 



20041 ). we began a small radial velocity monitoring program in 2004 with goal of finding addi- 
tional planets in already known planetary systems. Our methodology was to take advantage 
of the Hobby-Eberly Telescope's (HET) queue scheduled operations to obtain high cadence 
radial velocities of selected systems. In addition, we chose to obtain multiple radial velocity 
measurements per observational epoch to reduce the impact of both statistical and stellar 
noise on the search for previously hidden planets. This intensive technique limits the number 
of possible targets, but is a natural complement to the larger planet search programs. 

In this paper we report the first result from this project, the detection of a third planet 
around HD 74156. In §2 we give the properties and background information on the star. In 
§3 we describe our radial velocity measurements with the HET. We present our analysis that 
yields the detection of the new planet in §4. We conclude with a discussion of our result in 
^5. 



2. HD 74156 PROPERTIES 



According to the Hipparcos catolog flESAl 119971 ). HD 74156 (HIP 42723) is a GO star 
about 65 pc from the Sun {i^abs = 15.49±1.10 mas) with visual magnitude V = 7.61 and color 
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B — V = 0.59. iNaef et al.l (12004 ) previously reported the discovery of two planetary mass 
companions to the star based on its radial velocity variations. The planets had estimated 
minimum masses Mbsini = 1.9 Mj^p and M^sini = 6.2 Mj^p assuming the primary's mass 
M4 = 1.27 Mq. The lower mass planet is in a short period orbit {Pi, = 51.6 day), while the 
higher mass planet has a longer orbital period {Pc = 2025 day). Both planets' orbits were 
reported as highly eccentric (e = 0.64 and 0.58). 



Valenti fc Fischen (120051 ) de termined the i ron ab undance [Fe/H 



+0.13 for HD 74156. 

This result agrees with that from lSantos et al.l (j2003l ). who found [Fe/H] = +0.15. Therefore, 



it is metal ric h compared to the solar neighborhood, but average among host stars to high 
mass planets (IFischer &: Valenti II2005I ). 



Naef et al.l (120041 ) suggested that HD 74156 m ight be slightly evolv ed based on its large 



estimated luminosity relative to its spectral type. iTakeda et al.l (120071 ) report an age of 3.7 
± 0.4 Gyr based on a comparison of the star's properties with theoretical models. This age 
estimate implies that HD 74156 is sti ll on the m a in seq uence, although ages for field stars 
are notoriously difficult to determine. iNaef et al.l (120041 ) reported seeing no emission in the 
Ca II H line core, which is consistent with the star having at least an intermediate age. 



Takeda et all (120071 ) also proposed a mass Ma = 1.24 ± 0.04 Mq for HD 74156. This 
is consistent with the mas s adop ted in the discovery paper {Ma = 1.27 M©), which was 



estimated by ISantos et al.l (120031 ). We elected to adopt the former value in this paper for 



consis tency with our previous work. We included the uncertainty suggested by lTakeda et al. 
(120071 ) in the calculation of the planet minimum masses given in §4. 



RADIAL VELOCITY MEASUREMENTS 



We made high precisior i radial velocity m easurements of HD 74156 using the iodine 



absorption cell method (e.g. iButler et al.lll996l ). T he specific details of our implementation 



are the same as for our previous work described in iBean et al.l (|2007l ). 



Observations using the HET to feed the High Resolution Spectrograph (HRS. lTulllll998l ) 
were carried out by on-site telescope operators during 82 nights between UT dates 2004 
December 3 and 2007 May 11. The HRS was used in the resolution R = 60,000 mode 
with a 316 line mm~^ echelle grating. The cross dispersion grating was positioned so that 
the central wavelength of the order that fell in the break between the two CCD chips was 
5936 A. A temperature controlled cell containing molecular iodine gas (I2) was inserted in 
front of the spectrograph slit entrance during all exposures to imprint lines that provided 
a contemporaneous wavelength scale and instrumental profile fiducial for the high precision 
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radial velocity measurements. The exposure times were nominally 150 s, but varied up to 
three times that occasionally to account for increased seeing and/or cloud cover. Three 
sequential exposures were obtained each night and a total of 242 spectra of HD 74156 for 
radial velocity measurements were collected. 

HD 74156 was also observed once on 2006 December 14 without the iodine cell and with 
the same instrument setup, but in the R = 120,000 mode. The exposure time was 650 s. The 
spectrum from this observation served as a template for the radial velocity measurements. 

CCD reduction and optimal sp ectral extraction were carr ied out for all the individual 
spectra using the REDUCE package (jPiskunov &: Valenti II2002I ) and nightly calibration data. 
Relative radial velocities were measu r ed fro m the spectra using the modeling algorithm that 
is de scribed in detail by iBean et al.l (120071 ). We used the Gaussfit program (jjefferys et al. 



19881 ) to reduce each of the 82 multiple observation sets to a single velocity. The velocities we 
measured are relative to an arbitrary zero point, and we determined an offset simultaneously 
with our orbit analysis to adjust them to the system barycenter (see §4). The final velocities 
with this offset subtracted are given in Table [TJ They have a median uncertainty of 2.6 m 



4. ORBIT ANALYSIS 



We combined our measured radial velocities for HD 74156 with those from NaefetaL 



(12004 "ELODIE" and "CORALIE" samples) to create a dataset that spans 9.33 years. No 
other high precision radial velocities for HD 74156 have been published. The timespan, 
number of data points, and RMS residuals to the final best fit model for each of the three 
individual velocity samples are given in Table [2l The ELODIE and CORALIE measurements 
overlap, but there is a gap of 1.19 years between the final CORALIE measurement and our 
first measurement with the HET. 

We fit the total radial velocity dataset with a double Keplerian model to determine 
the two previously known planets' period, P, velocity semiamplitude, K, eccentricity, e, 
longitude of periastron, u, and time of periastron, Tp. To account for the heterogeneous 
nature of the dataset, we also simultaneously determined offset parameters for each individual 
velocity sample to adjust them to be relative to the sy stem barycenter. We performed the 
fitting using the Gaussfit program (jJefferys et al.lll988l ) with both robust and least squares 
("reduced chi squared" metric, xl) estimation to determine the parameter values that gave 
best match between our model and the measured data. For our robust fitting we used the 
"fair" metric described by iRey Jl983h with an adopted asymptotic relative efficiency (ARE) 



- 5 - 



of 0.92. We refer to the models found with the two methods as "Robust" and "Least Squares" 
hereafter. In the case of the two planet model, the parameters determined using the two 
methods were consistent. 

To searc h for evidence of additional planets in the HD 74156 system we calculated a 
periodogram (jPress et al.lll992l ) of the velocity residuals from the two planet model fit. The 
result is shown in Figure [H The highest peak in the periodogram is at a period around 349 
days and has a form al false alarm pro bability (FAP) of 0.0014%. We used the bootstrap 
method described by lEndl et al.l (120021 ) to further investigate the likelihood that this signal 
was real. We calculated a periodogram and noted the maximum power for each of 10,000 
simulated data sets. The data sets were generated by randomly scrambling the two planet 
model residuals while maintaining the sampling of the observations. We did not find a 
maximum power in any of the trial periodograms as high or higher than the peak at 349 
days in the periodogram of the real two planet model residuals (19.98). This implies a FAP 
for the 349 day signal of less than 0.01%, which is consistent with the formal estimate. We 
took this as evidence for a real periodic signal in the radial velocities not accounted for by 
the two previously known planets. 

We explored the possibility that the detected periodic signal is attributable to a third 
planet in the system by fitting a triple Keplerian model to the radial velocity dataset. We 
again used both the robust and least squares estimation methods. The quality of the fits 
improved significantly with the addition of the third component. However, we found that 
the RMS of the HET velocity residuals (6.0 m s~^) remained more than twice the median 
of the measurement uncertainties (2.6 m s^^). This discrepancy could be due to using an 
incomplete model of the system (i.e. there are more than 3 planets), intrinsic variations of 
the star's photosphere that mimic real radial velocity changes, improper weighting of the 
ELODIE a nd CORAL IE d ata, and/or e r rors in our velocity measurement method. Using the 
method of IWrightl (120051 ). iButler et al.l (120061 ) estimated the radial velocity "jitter" for HD 
74156 to be 4.0 m s~^. Adding this value in quadrature with the HET velocity uncertainties 
yields errors that are slightly above what would be expected from the fit residuals and the fit 
xl = 0.88. This is not surprising because our method of taking multiple measurements over 
10 - 25 minutes and combining them into one measurement reduces the noise arising from 
short-term stellar variability. We ultimately decided to increase the HET uncertainties by a 
factor of two in order to account for potential errors in the data and to avoid overweighting 
them in the orbital fitting, but we cannot discount the possibility that there are additional 
planets in the system causing the higher than expected dispersion. Therefore, we have elected 
to publish our original estimate of the HET velocity uncertainties in Tabled] to avoid biasing 
future investigations. 
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We repeated both the two and three planet orbit fitting using the revised uncertainties 
for the HET data. The orbital parameters and associated uncertainties that we determined 
for the three planet models and minimum masses calculated from these parameters are given 
in Table O We found that the xl of the Least Squares model improved from 2.3 to 1.4 
by introducing the third planet into the model. The RMS of the HET velocity residuals 
from the fits improved from 8.5 to 6.0 m s~^ for both methods. The RMS for the ELODIE 
residuals went from 12.7 to 10.8 m s~^ in the Least Squares model and from 12.8 to 9.5 m 
s~^ in the Robust model. However the RMS of the CORALIE residuals degraded from 10.3 
to 10.5 m s~^ in the Least Squares model and from 10.2 to 11.1 m s~^ in the Robust model. 
The three planet model does yield an improved fit to the combined ELODIE and CORALIE 
dataset, with the RMS of the residuals from the Robust model improving from 11.7 to 10.3 
m s~^. 

The HET velocities and the Robust fit are shown in Figure [2l The radial velocity 
data is shown phased to each component's period and with the other two components' 
orbits subtracted assuming the Robust model in Figure [31 The Robust and Least Squares 
models for the two previously known planets are very similar. However, the period, time of 
periastron, and eccentricity for the candidate third planet determined using the two different 
methods are inconsistent. The most notable difference between the models for this planet is 
in the derived eccentricities, with the Robust model having e = 0.25 and the Least Squares 
model having e = 0.55. For comparison, the new planet's phased orbit for the Least Squares 
model is shown in Figure |H 

The least squares method of fitting involves weighting the data points according to 
their input variances and is based on the assumption that the data are normally distributed. 
The purpose of using robust estimation instead is to reduce the influence of non-normal 
outliers in the data. Rather than attempting to identify outliers by hand, robust metrics 
place outlier identification and their subsequent reweighting on a rigorous footing. Data 
points with fit residuals larger than a certain threshold are iteratively down-weighted, rather 
than completely ignored, during the fitting process. However, the identification of some 
data points as outliers does not mean they are necessarily bad. It only means that they 
don't agree with adopted physical and noise models and thus should not carry their normal 
weight when assessing the fit quality to avoid avoid biasing the final solution. For example, 
a small signal due to a fourth planet in the system could cause some data points to seem 
highly discrepant. In this case the data points would be correct, but would adversely affect 
the fitting of a model that only includes three componenets if they were not identified and 
down-weighted. 

The similarity of the derived orbital parameters for the two previously known planets 
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when using the robust and least squares methods is due to the fact that even the worst 
possible outliers in the radial velocity data are small relative to the amplitudes of the plan- 
ets' velocity signals. However, the size of the outliers is the same or even larger than the 
amplitude of the velocity signal from the possible third planet. Therefore, whether potential 
outliers are identified and treated differently significantly influences the orbital solution for 
this component. 

The Robust model for HD 74156 d is based on the estimation that some of the HET 
data points are actually contaminating outliers. This isn't surprising given the variety of 
physical processes that can influence a measured stellar radial velocity at the 2 - 3 m s~^ 
level. The down- weighting of the identified outliers leads to a less eccentric orbit as the fit to 
highly deviant points with small input uncertainties is relaxed. The Least Squares model is 
the best fit to the data assuming that the input errors are scaled correctly and the data are 
normally distributed. The weight given to all the HET data points is based on their input 
uncertainties, which are twice those given in Table [H The result is a more eccentric orbit 
for the new planet because the solution is strongly influenced by the potential outliers. We 
favor the Robust model over the Least Squares model, but we chose to present and discuss 
both in order to illustrate the uncertainty in derived orbital parameters for planets when the 
signal is near the detection threshold in the data. 

A periodogram of the three planet Robust model residuals is shown in Figure [51 The 
periodogram of the residuals to the Least Squares model is essentially the same. The addition 
of the third component to the model removes the periodicity detected at 349 days. No 
periodicity is detected in the three planet model residuals with a FAP less than 22%. This 
indicates that a summation of three Keplerian orbits adequately accounts for the detected 
periodic signals in the radial velocities. 

We note that the periodogram of the two planet model residuals in Figured] also exhibits 
moderately significant power around 175 days. This is roughly half the value of the period we 
explored and ultimately determined for the potential third planet. Also, the 175 day spike 
in the two planet model residual periodogram has a smaller power (FAP = 0.09% from the 
bootstrap simulation) than the spike at 349 days and is not present in the periodogram of the 
three planet model residuals. We therefore suspect that the smaller, moderately significant 
spikes around 175 days in the periodogram of the two planet model residuals were aliases. 
Nevertheless, we attempted to fit a three planet model to the radial velocity data where the 
third component had a period around 175 days. We were not able to determine a physically 
realistic set of parameters for the third component with this approach and various trials 
with fixed parameters yielded much worse fits than the three planet fits where the third 
component had a period close to 349 days (A^^ > 80 for the least squares approach). 
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The significant improvement in fit quality for botli tlie Robust and Least Squares models 
and comparison of the two planet and three planet model residual periodograms lead us to 
propose that there is a third planet in the HD 74156 system with orbital parameters similar 
to those given in Table O Beyond the detection of this new planet, our data and analysis has 
also yielded improved orbital parameters for HD 74156 c. We find for t his planet a perio d 



22% longer and eccentricity reduced by 0.13 from the values determined by lNaef et al.l (120041 ) . 
This is not an unusual result when determining revised orbital parameters for long period 
planets with the addition of significantly more radial velocity data. 

From our determined Robust model parameters we estimate Mbsini = 1.88 ± 0.03, 
Mrsini = 8.03 ± .12, and Mrfsini = 0.40 ± 0.02 Mjup assuming Ma = 1.24 ± 0.04 Mq 



( iTakeda et al.ll2007l ). The derived velocity semiamplitude for the new planet (~ 11.5 m s^^) 
is roughly the same magnitude as the ELODIE and CORALIE residuals, while it is roughly 
twice that of the HET residuals. Therefore, the HET velocities are the most sensitive to the 
third planet. Its period is also close to a year (but significantly different enough so that it 
cannot be due an alias or error in the data) and the 2.5 years of HET observations do not 
give full phase coverage. More high precision observations over the next few years will be 
needed to refine the orbital parameters for this planet. 

The eccentricities we have determined by fitting the radial velocity data indicate that 
the two previously known planets are in highly eccentric orbits, while the new planet is in 
at least a moderately eccentric orbit. This raises the question of whether the derived orbital 
parameters represent a stable three planet system. We made a preliminary check of this 
using the Runge Kutta numerical integrator in the Systemic Consol^ We integrated the 
planet positions forward in time for 1000 years assuming their true masses were equal to the 
minimum masses and that their orbits were in the same plane. We used a 0.05 day time 
step for the calculations. The 1000 year integration represents about 7100 full orbits of the 
innermost planet. We found that the semi-major axes of all three planets were not predicted 
to change from their initial values during the entirety of both integrations. Therefore, both 
sets of orbital parameters that we have determined do not describe highly unstable systems. 



Prior to our discovery, [Raymond &: Barnes! (120051 ) had determined that an additional 



Saturn-mass planet in the HD 74156 system having an orbital semimajor axis a = 0.9-1.4 
AU, and e < 0.15 would very likely be dynamically stable for 100 million years. They had 
even predicted that because such a plan et could exist, then it mu st exist based on their 



"Packed Planetary Systems" hypothesis (IBarnes fc Raymond! |2004| ). The candidate planet 



we have identified has M sini = 1.3 Msat and a = 1.0 AU, which is very similar to the prop- 



Available at 



http : // www . oklo . org 



We used the May 30, 2007 version of the code. 
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erties of their stable test planet. However, the simulations presented by lRaymond &: Barnes 
( 120051 ) were based on orbital parameters for the previously known planets that are some- 
what different than we have determined based on the new radial velocity data and we find 
a slightly higher eccentricity for the new third planet than they considered. Nevertheless, 
the results of their robust simulations support our proposal that the third planet exists and 
that the orbital parameters we have derived for it are plausible. 



5. DISCUSSION 



With our detection of a third planet around HD 74156, it becomes the eighth normal 
star to host three or more planets. Currently, five other stars are known to host three planets 
(GJ 876, V And, GJ 581, HD 69830, and HD 37124) and two are known to host four planets 
(p^ Cnc and /i Arae). The most recently determined orbital eccentricities and periods for 
these planets are tabulated along with the source of the data in Table HI These same data 
are plotted in Figure O The planets in systems containing three or more planets have a 
median orbital eccentricity of 0.13, which is significantly lower than the median of 0.25 for 
all known exoplanets. This could be due to the increased likelihood of dynamical instabilities 
to develop in higher ord er planetary systerns whe re one or more component has a significantly 
elliptical orbit. Indeed, IChatterjee et al.l (120071 ) have proposed that dynamical instabilities 
leading to planet ejection can quickly arise in systems initially containing three gas giants 
without requiring special orbital configurations or additional evolutionary mechanisms. It 
might then be expected that the three planet systems which survive for the timescale of 
the typical known exoplanet host star's age (i.e. those that are detectable with current 
techniques) would be those with planets having low orbital eccentricities. 

This hypothesis seems to be supported by the sample of previously known systems con- 
taining three or more planets. However, HD 74156 b and c are in orbits that are much 
more eccentric than any of the components in these systems. Is HD 74156 a rare exam- 
ple of a system that survived a period of dynamical instability with three gas giants in 
highly eccentric orbits, or does its existence point to other mechanisms as playing an im- 
portant role in the evolution of planetary systems? Studies of the dynamical characteris- 



tics of HD 74156 under the assumi 


otion of the previous two p 


anet model have been car- 


ried out by \ 


Vaeasawa et al. (2 


003). 


Barnes & Ravmond (2004), 


Ravmond & Barnes (2005), 


Ravmond et al. (20061). Adams & Lauehlin (2006). Barnes & Greenbere (2006). Libert & Henrad 


(2006a), and 


Libert & Henrad 


(2006bh. Clearlv. new dvnamical studies should be undertaken 



now that we have uncovered the third planet in order to understand the possible origins of 
this system's unique configuration. 



-10- 



We would like to thank Eugenio Rivera, Stefano Meschiari, Aaron Wolf, and Greg 
Laughlin for providing the Systemic Console to the community. We also thank an anony- 
mous referee for a careful reading of the manuscript and constructive comments which im- 
proved this paper. Wc acknowledge support from NASA GO-09407, GO-09969, GO-09971, 
GO-10103, GO-10610, and GO-10989 from the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555; and from JPL 1227563 (S'/MMASSIF Key Project, Todd Henry, P.I.), 
administered by the Jet Propulsion Laboratory. 



- 11 - 



REFERENCES 

Adams, F., & Laughlin, G. 2006, ApJ, 649, 992 
Barnes, R., & Greenberg, R. 2006, ApJ, 652, L53 
Barnes, R., & Raymond, S. N. 2004, ApJ, 617, 569 

Bean, J. L., McArthur, B. E., Benedict, G. F., Harrison, T. E., Bizyaev, D., Nelan, E., &; 
Smith, V. V. 2007, AJ, 134, 749 

Benedict, G. F., et al. 2002, ApJ, 581, 115 

Benedict, G. F., et al. 2006, AJ, 132, 2206 

Bouchy, F., et al. 2005, A&A, 444, L15 

Butler, R. R, Marcy, G. W., Williams, E., McCarthy, C., Dosanjh, R, & Vogt, S. S. 1996, 
PASP, 108, 550 

Butler, R. P., et al. 2006, ApJ, 646, 505 

Castellano, T. 2000, PASP, 112, 821 

Charboneau, D., Brown, T., M., Latham, D. W., & Mayor, M. 2000, ApJ, 529, L45 



Chatterjee, S., Ford, E. B., & Rasio, ApJsubmitted, astro-ph/0703166 
Chiang, E., & Murray, N. 2002, ApJ, 576, 473 

Cowan, N. B., Agol, E., & Charbonneau, D. 2007, MNRAS, 379, 641 

Deming, D., Seager, S., Richardson, L. J., & Harrington, J. 2005, Nature, 434, 740 

Endl, M., Kiirster, M., Els, S, Hatzes, A. P., Cochran, W. D., Dennerl, & K. Dobereiner, S. 
2002, A&A, 392, 671 

ESA 1997, The Hzpparcos and Tycho Catalogues (ESA SP-1200; Noordwijk: ESA) 

Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 

Ford, E. B., Lystad, V., & Rasio, F. A. 2005, Nature, 434, 873 

Gillon, M., et al. 2007, A&Aaccepted, astro-ph/0705.2219 

Harrington, J., et al. 2005, Science, 314, 623 



- 12 - 



Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, ApJ, 529, L41 
Janson, M., et al. 2007, AJ, 133, 244 

JefFerys, W. H., Fitzpatrick, M. J., & McArthur, B. 1988, Celestial Mechanics 41, 39 
Knutson, H. A., et al. 2007, Nature, 447, 183 
Laughlin, G., & Adams, F. C. 1999, ApJ, 526, 881 
Libert, A.-S., & Henrad, J. 2006a, A&A, 461, 759 
Libert, A.-S., & Henrad, J. 2006b, Icarus, 183, 186 

Lopez-Morales, M., Morrell, N. I., Butler, R. R, & Seager, S. 2006, RASP, 118, 1506 
Lovis, C., et al. 2006, Nature, 441, 305 
McArthur, B. E., et al. 2004, ApJ, 614, L81 

Naef, D., Mayor, M., Beuzit, J. L., Perrier, C, Queloz, D., Sivan, J. P., & Udry, S. 2004, 
A&A, 414, 351 

Nagasawa, M., Lin, D. N. C., & Ida, S. 2003, ApJ, 586, 1374 

Pepe, F., et al. 2007, A&A, 462, 769 

Piskunov, N. E., & Valenti, J. A. 2002, A&A, 385, 1095 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes: 
The Art of Scientific Computing (2nd Edition; Cambridge, U.K.: Cambridge Univer- 
sity Press) 

Raymond, S. N., & Barnes, R. 2005, ApJ, 619, 549 

Raymond, S. N., Barnes, R., & Kaib, N. A. 2006, ApJ, 644, 1223 

Rey, W. J. J. 1983, Introduction to Robust and Quasi-Robust Statistical Methods (New 
York: Springer- Verlag) 

Richardson, L. J., Deming, D., Horning, K., Seager, S., & Harrington, J. 2007, Nature, 445, 
892 

Rivera, E. J., et al. 2005, 634, 625 
Sato, B., et al. 2005, ApJ, 633, 465 



Santos, N. C, Israelian, G., Mayor, M., Rebolo, R., & Udry, S. 2003, A&A, 398, 363 
Shankland, P. D., et al. 2006, ApJ, 653, 700 

Takeda, G., Ford, E. B., Sills, A., Rasio, F. A., Fischer, D. A., & Valenti, J. A. 2007, ApJS, 
168, 297 

TuU, R. G. 1998, Proc. SPIE, 3355, 387 

Udry, S. et al. 2007, A&A, 469, L43 

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 

Wright, J. T. 2005, PASP, 117, 657 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-14- 



Table 1. HEX Radial Velocities for HD 74156 



HJD - 2450000.0 RV (m s"!) 



3342, 


,8969 


-47, 


,28 


± 


3, 


,84 


3347, 


,0019 


-41, 


,26 


± 


2, 


,.58 


3355, 


,8330 


-34, 


,26 


± 


2, 


,89 


3357, 


,8444 


-35, 


,63 


± 


2, 


,91 


3359, 


,8504 


-28, 


,42 


± 


4, 


,95 


3360, 


,9748 


-37, 


,04 


± 


3, 


,19 


3364, 


,9763 


-41, 


,36 


± 


3, 


,23 


3365, 


,8259 


-40, 


,73 


± 


2, 


,89 


3367, 


,8213 


-64, 


,14 


± 


2, 


,93 


3383, 


,9213 


-64, 


,01 


± 


3, 


,37 


3390, 


,7526 


-27, 


,69 


± 


2, 


,82 


3448, 


,7387 


33, 


,62 


± 


1, 


,59 


3451, 


,7306 


27, 


,48 


± 


2, 


,41 


3476, 


,6500 


-118, 


,87 


± 


2, 


,57 


3480, 


,6372 


-127, 


,55 


± 


5, 


,68 


3481, 


,6318 


-82, 


,68 


± 


2, 


,08 


3482, 


,6317 


-57, 


,03 


± 


2, 


,13 


3664, 


,9959 


141, 


,47 


± 


1, 


,96 


3675, 


,9704 


127, 


,44 


± 


2, 


,24 


3676, 


,9844 


123, 


,13 


± 


3, 


,69 


3682, 


,9515 


20, 


,69 


± 


2, 


,78 


3687, 


,9321 


15, 


,18 


± 


2, 


,50 


3689, 


,9285 


61, 


,28 


± 


2, 


,54 


3691, 


,9156 


92, 


,73 


± 


2, 


,98 


3697, 


,9126 


123, 


,87 


± 


2, 


,03 


3703, 


,8858 


147, 


,52 


± 


4, 


.60 


3708, 


,8817 


146, 


,26 


± 


3, 


.24 


3710, 


,8789 


150, 


,70 


± 


2, 


,58 


3718, 


,0141 


136, 


,36 


± 


4, 


,98 


3724, 


,8303 


129, 


,03 


± 


3, 


,12 


3728, 


,8252 


121, 


,82 


± 


2, 


.94 


3731, 


,9672 


83, 


,91 


± 


3, 


.86 


3733, 


,8020 


39, 


,79 


± 


3, 


.20 


3734, 


,8087 


-3 


.84 


± 


2, 


,64 


3736, 


,9455 


-86, 


,32 


± 


3, 


,54 


3741, 


,7833 


60, 


,98 


± 


2, 


,69 


3742, 


,7830 


73, 


,14 


± 


2, 


,33 


3743, 


,7904 


83, 


,45 


± 


2, 


.72 


3748, 


,7726 


128, 


,14 


± 


3, 


.04 


3751, 


,7677 


139, 


,71 


± 


3, 


.00 


3753, 


,7634 


150, 


,89 


± 


2, 


,61 


3754, 


,7496 


141, 


,44 


± 


2, 


,43 


3756, 


,7485 


141, 


,78 


± 


3, 


,11 


3764, 


,7353 


140, 


,77 


± 


2, 


,69 


3832, 


,6754 


100, 


,92 


± 


2, 


.17 


3833, 


,6962 


94, 


,87 


± 


2, 


.55 
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Table 1 — Continued 



HJD - 2450000.0 RV (m s"!) 



3834, 


,6755 


85, 


,38 


± 


1.92 


3835, 


,6671 


63, 


,55 


± 


2.01 


3838, 


,6638 


-35, 


.39 


± 


1.96 


3841, 


,6443 


-46, 


.03 


± 


4.92 


3845, 


.6318 


72, 


,93 


± 


2.72 


3846, 


,6533 


77, 


,92 




2.69 


4029, 


,9895 


134, 


,42 




2.58 


4035, 


,9910 


116, 


.77 




3.07 


4038, 


,9782 


102, 


.21 


± 


1.65 


4039, 


,9707 


102, 


.13 


± 


2.46 


4040, 


.9590 


79, 


.20 




2.02 


4043, 


.9637 


11, 


.48 


± 


1.25 


4044, 


.9538 


-18, 


.63 


± 


1.81 


4045, 


.9556 


-72, 


.29 


± 


2.05 


4050, 


.9539 


42, 


.33 


± 


2.30 


4051, 


.9467 


52, 


.82 


± 


2.67 


4052, 


.9393 


67, 


.80 


± 


1.72 


4073, 


.8806 


125, 


.65 


± 


1.92 


4079, 


.8634 


130, 


.08 


± 


2.64 


4087, 


.8437 


95, 


.94 


± 


1.78 


4106, 


.7835 


63, 


.79 


± 


2.55 


4109, 


.7887 


93, 


.54 


± 


2.53 


4110, 


.7994 


106, 


.19 


± 


2.29 


4129, 


.8706 


105, 


.98 


± 


2.14 


4130, 


.7409 


103, 


.57 


± 


2.36 


4133, 


.8470 


88, 


.55 


± 


2.69 


4134, 


.7264 


92, 


.38 


± 


2.58 


4135, 


.8673 


103, 


.20 


± 


2.85 


4136, 


.8409 


90, 


.93 


± 


3.26 


4148, 


.6771 


-81, 


.12 


± 


3.41 


4156, 


.6604 


48, 


.73 


± 


3.71 


4159, 


.7804 


66, 


.51 




2.32 


4166, 


.7627 


89, 


.64 


± 


2.10 


4167, 


.7585 


85, 


.27 


± 


2.00 


4211, 


.6299 


67, 


.62 


± 


2.47 


4231, 


.6001 


104, 


.84 


± 


1.78 
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Table 2. The Radial Velocity Samples 



Sample 


Time Span 


N 


RMS (m 8-1)=- 


ELODIE 


1998.02 


- 2003.37 


48 


9.5 


CORALIE 


2001.03 


- 2003.74 


38 


11.1 


HET 


2004.92 


- 2007.36 


82 


6.0 



''RMS of the residuals from the Robust three planet 
model. 
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Table 3. Derived Parameters for the HD 74156 Planets 



Parameter HD 74156 b HD 74156 c HD 74156 d 
Robust solution 

P (days) 51.65 ± 0.01 2476.7 ± 8.7 346.6 ± 3.6 

K(ms-i) 115.1 ± 1.5 115.5 ± 2.3 10.5 ± 1.2 

To (HJD- 2450000.0) 1980.8 ± 0.1 952.2 ± 13.8 678.2 ± 44.2 

e 0.64 ± 0.01 0.43 ± 0.01 0.25 ± 0.11 

uj (deg) 175.8 ± 0.8 261.3 ± 2.0 166.5 ± 27.4 

Msini {Mj„p)^ 1.88 ± 0.03 8.03 ± 0.12 0.40 ± 0.02 

Least Squares solution 

P (days) 51.65 ± 0.01 2481.8 ± 8.5 339.6 ± 2.6 

if(ms-i) 114.9 ± 1.8 118.1 ± 2.8 13.5 ± 2.0 

To (HJD-2450000.0) 1980.8 ± 0.1 927.6 ± 14.3 757.5 ± 25.7 

e 0.64 ± 0.01 0.44 ± 0.01 0.55 ± 0.07 

uj (deg) 175.8 ± 0.8 257.6 ± 2.3 188.1 ± 11.1 

Msini{Mj„p)^ 1.87 ± 0.03 8.19 ± 0.13 0.45 ± 0.04 



Assuming Ah = 1.24 ± 0.04 Mq l lTakeda et al.ll2007h . 
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Table 4. Orbital Parameters for Planets in Systems with Three or More Planets 



Planet P (days) e Source 



HD 74156 b 


51.65 





64 


1 


HD 74156 c 


2476.7 





43 


1 


HD 74156 d 


346.6 





25 


1 


GJ 876 b 


60.83 





03 


2 


GJ 876 c 


30.46 





26 


2 


GJ 876 d 


1.94 





00 


2 


pi One b 


14.65 





01 


3 


pi One c 


44.36 





07 


3 


pi One d 


5552 





09 


3 


pi One e 


2.80 





09 


3 


V And b 


4.62 





02 


3 


V And c 


241.2 





26 


3 


V And d 


1290.1 





26 


3 


GJ 581 b 


5.37 





02 


4 


GJ 581 c 


12.93 





16 


4 


GJ 581 d 


83.6 





20 


4 


HD 69830 b 


8.67 





10 


5 


HD 69830 c 


31.56 





13 


5 


HD 69830 d 


197.0 





07 


5 


p Arae b 


643.25 





13 


6 


p Arae c 


9.64 





17 


6 


p Arae d 


310.55 





07 


6 


p Arae e 


4205.8 





10 


6 


HD 37124 b 


154.46 





06 


3 


HD 37124 c 


2295 





20 


3 


HD 37124 d 


843.6 





14 


3 



R eferences. — (1 ) This P aper, (2) 
iRlvcra ct all ||2005|). (3) iButler et al, (200§,), 



4) lUdrv et all li2007h. [ 5) iLovis et al.l 
2006h . f6') |Pepe et all ll2"o07h 
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Fig. 1. — Periodogram of the velocity residuals from the two planet Robust model. The 
spike around 349 days indicates a remaining periodic signal that is not accounted for by the 
previously known planets. The false alarm probabihty of the peak is 0.0014%. 
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Fig. 2. — Top: HET radial velocities (circles) and the Robust three planet model to all the 
data (line). Bottom: Residuals from the fit (circles). The error bars are omitted for clarity. 
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Fig. 3. — Representation of the Robust three planet model. The radial velocity data (circles) 
are phased to each component's period with the other two components' orbits subtracted 
and the component's model velocities (lines). The error bars are omitted for clarity. 
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Fig. 4. — Representation of the Least Squares orbit for the third planet (hne) with the other 
two components' orbits subtracted from the data (circles). Note that the small differences 
in the parameters for the two other planets do yield slightly different residuals when their 
orbits are subtracted from the data. The error bars are omitted for clarity. 
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Fig. 5. — Periodogram of the residuals to the three planet Robust model. No periodicity is 
detected with false alarm probability less than 22%. 
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Fig. 6. — Orbital eccentricity as a function of period for the planets in the eight systems 
containing three or more planets (filled circles, see Table H] for the data and sources). The 
median value for this sample is 0.13 (indicated by the solid line), while the median orbital 
eccentricity of all the known exoplanets is 0.25 (indicated by the dashed line). The planets 
in the HD 74156 system (open circles) are notably eccentric relative to the planets in the 
systems with three or more planets. Note that the orbital eccentricity of HD 74156 d (P = 
347 days) is poorly constrained with the current data. 



